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Abstract 

The effect of the photon emission (bremsstrahlung) in the cross section of the process 
of direct production of the Higgs boson in the future high luminosity electron and 
muon colliders is calculated. It was found that cross section at the top of the Higgs 
boson resonance peak is reduced by factor 0.347 for the electron collider and 0.548 
for the muon collider. Machine spread of the centre of the mass energy of 4.2MeV 
(equal to the Higgs width) would reduce peak cross section further, by factor 0.170 
and 0.256 (QED and energy spread) for electron and muon beams respectively. 
Possible uncertainties in the resummed QED calculations are discussed. Numerical 
results for the lineshape cross section including QED and many values of the machine 
energy spread are provided. 
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1 Introduction 


Observation of the Higgs particle at CERN Si has opened the new era of precise study of 
its properties. Once its mass is known, one may reconsider a possibility of its direct obser¬ 
vation in the lepton annihilation e“e+ H or fjT —)■ H, at future electron colliders or 
dedicated muon collider. Due to extremely small coupling of the Higgs boson to electron, 
it seems at hrst sight that its direct production in the electron collider is just hopeless. 
However, in the Future Circular Collider with beams (FCCee) considered at CERN 
featuring very high luminosity, this process would in principle be observable, provided one 
could eliminate copious background processes. On the other hand, stronger coupling of 
Higgs to muons puts dedicated muon collider into a dehnite advantage, provided decent 
luminosity and small beam energy spread are achieved. In either case, machine energy 
spread and an additional smearing of the beam energy due to QED bremsstrahlung are 
major points in the feasibility studies of these projects. This is why the present study was 
undertaken. The influence of the beam energy spread on observability of the direct Higgs 
observation in FCCee collider was already discussed at the 8th FCC-ee Physics Work¬ 
shop [^. In this work we shall concentrate mainly on calculating effects due to initial 
state radiation of multiple photons. This will be done exploiting past experience in calcu¬ 
lating very precisely similar QED effect for the Z boson productions at LEP experiments, 
for instance in ref. [^. Similar analysis of the intial state QED corrections, taking into 
accound the machine energy spread, for —?■ H process can be also founc0 in refs. 

and (^. 

The paper is organized as follows: After defining the Born cross section for the Higgs 
process production in the lepton annihilation, we will discuss the effect of the initial state 
radiation (ISR) corrections. We shall discuss theoretical uncertainties in the evaluation 
of this QED effect, presenting numerical results for three different QED ISR formulas, 
of the varying sophistication level. The effect of the machine energy spread in the cross 
section will be added in the discussion, presenting numerical results for several values of 
such a spread. Finally the analogous effects in planned muon colliders will be discussed 
and numerical results will be presented. 


2 QED initial state radiation formulas 


Born cross section for the Higgs particle production in collider is given by the (rela¬ 
tivistic [^) Breit-Wigner (B-W) formula 


0-5(5) 


{s - Mlf + 


( 2 . 1 ) 


where Mh = 125.7 GeV and Fh = 4.2MeR according to refs. [9|10 
~e~ —)■ is Bee = 5.3 ■ 10“® 


The branching ratio 

of —)■ is Bee = 5.3 ■ 10“^. The electron and muon branching ratios are related by 

the factor m't/ml, thus the above result can be obtained from B^^ = 2.19 ■ 10“^, see also 


^We thank M. Greco and S. Dittmaier for bringing these works to our attention. 
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refs. [ 9 , 10 . There are variants of the B-W formulas with the s-dependent width, but for 


the narrow resonance like Higgs they differ negligibly from the above, see more discussion 
in the Appendix A. 

The initial state radiation correction to this process was calculated using formulas 
of ref. for Z boson productioij^ The entire initial-state 0( 
integrated cross section reads: 


a 


/prag 


formula of ref. 


aj(s) = J dv pi{v)aB{s{l - v)), 

0 

pi{v) = F( 7 ) {ds Ah{v)} , 


( 2 . 2 ) 


where 


I 


(aj 

ip) 


da 


1 + 1 + 

1 + i + ^ + Af 


Aniv) 


(-1 + i) +7 


0 

1 + i) 

l+3(l-t>)^ 


(2.3) 


ln(l 


and 


OYFS — 7 + - —o + ^ 


7 = 2f(ln;^-l), ^ (7) = 


(2.4) 


Here, a is the QED coupling constant, irtf. the electron mass and C is the Euler-Mascheroni 
constant. In case of the muon beams rrie in 7 is replaced by and i?ee by 

Zero spin nature of the Higgs boson instead of spin one of Z counts negligibly in the 
QED ISR effects, simply because deformation of the resonance curve is mainly due to soft 
photons. The constant A which is responsible for the above spin difference is of order 
^ ~ 1/400 without any logarithmic enhancement. It is influencing mainly an overall 
normalization - hence at the precision level we are aiming at, it can be safely set to zerc]^ 
On the other hand soft photon exponentiation/resummation in eq. (2.2) is critical and 
mandatory. The formula of eq. (2.2) comes from standard diagrammatic perturbative 


QED calculations including Yennie-Frautschi-Suura (YES) exponentiation, see ref. 
and was originally introduced for the purpose of the algebraic validation of the Monte 
Carlo program YFS2 of ref. 11 . Later on it was discussed and used in many papers, see 
for instance refs. [5,12 -15 and the references therein. 


Three variants (for / = (a), (6) or (c)) of the ISR formula in eq. (2.2) correspond 
to the increasing sophistication (perturbative order) of the non-soft collinear radiative 
corrections. Change of the type of ISR formula will be used to estimate uncertainty due 
to unknown/neglected QED higher orders. 


^See eq. (202) therein. 

^ Constant A is also set to zero in ref. [^, while in ref. [^ vertex and real-soft contributions are 


provided, but non-logarithmic constant A is was explicitly obtained. 
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E 

(Jb 

(^(a) 

o-(fe) 

^(0 

O-R 

O-R 

O-R 

Mh 

1.6413 

0.5398 

0.5693 

0.5701 

0.328 

0.346 

0.347 

Mh + Eh 

0.3283 

0.1846 

0.195 

0.1949 

0.562 

0.594 

0.593 

Mh-Eh 

0.3283 

0.1073 

0.1134 

0.1133 

0.326 

0.345 

0.345 


Table 1: Born cross sections (in fb) of eq. (2.1) and three variants of the ISR-corrected cross 
sections of eq. (2.2) for three values of y/s = E = Mh = 125.7 GeV and E = Mh ± Eh in 
the electron-positron collider. The integration error is below 0.0005/6. 


3 Machine energy spread 


In real accelerator experiments the beam is not monoenergetic, i.e., centre-of-mass energy 
E = y/s has spread 6 around the centre value Eq = y/^ of the beam energy. The 
distribution of E is usually well approximated by the following Gaussian distribution 


G{E-Eq-5) = 


1 




(g--En) 


(3.1) 


In case of no QED effects, the Born cross section (2.1) gets simply convoluted with the 


energy spectrum of eq. (3.1): 


a 


B 


\E-6)= dE'aB{E')G{E'- E-6). 


(3.2) 


Once QED bremsstrahlung is switched on, the following double convolution provides 
realistic experimental cross section: 


^convi^E- 6) = j dE' ai{E')G{E' - E; 6) 

r /•! 


= / dE' 


/ dv , _ 

Iq 


^ Pi(v)<Tb(E'‘(1 - !>)). 


(3.3) 


for three variants, / = (a), (6), (c) of the radiative function (2.2) 


Because of rapid decrease of the Gaussian distribution for large arguments, the energy 
integration range will be restricted to E — 106 < E' < E + 106 without any loss of the 
calculation reliability. The numerical integrations in one and two dimensions require a 
little bit of care, because of strongly singular integrands. The adaptive integration library 
functions of ROOT library pT] were used. All results were also cross-checked using FOAM 


adaptive Monte-Garlo simulator/integrator of 17-19 


4 Numerical results for electron-positron colliders 


The Born cross sections and results from three variants of the ISR formula of eq. (2.2) 
are presented in Fig. The case and a^c) are almost indistinguishable. In addition, 

^ Integration errors were also taken from FOAM, as they are more reliable. 
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E 

ctb 

(1) 

(2) 

(3) 

w 

O-R 

M 

Mh 

1.6413 

0.5701 

0.2795 

0.1825 

0.347 

0.170 

0.111 

Mh + Eh 

0.3283 

0.1950 

0.2323 

0.1753 

0.594 

0.708 

0.534 

Mh-Eh 

0.3283 

0.1133 

0.1865 

0.155 

0.345 

0.568 

0.472 


Table 2: Cross section of the direct Higgs production process in the electro-positron collider 
including effect of ISR type (c) and the machine energy spread. No energy spread is in 
while (7(2) and (j( 3 ) include ISR and energy spread 6 = 4.2 MeV and 6 = 8.0 MeV according 
to eq. (3.3). Reference Born cross section cr^ is also included and the ratios with respect to 
Born are also provided. The integration error is below 0.0005/6. 


Tab. [T] presents the same results for the energy at the peak E = Mh and near the peak 
E = Mh ± Eh with the 4-digit precision. The QED uncertainty from unaccounted QED 
higher orders can be estimated from Tab. looking into differences between ISR type 
/ = (c) and I = ip). It is below 0.3% and is compatible with the estimate of neglected A 
in eq. (2.2) (it is also comparable to the numerical integration error). From now on we 
shall use ISR formula for I = (c) only. 

The introductory exercise on the machine energy spread is shown in Fig. where Higgs 
production cross section in the electron-positron collider is plotted without energy spread 
(Born) and for several values of the machine energy spread 6 = 4.2, 8,15, 30, lOOMeV. 
The QED ISR effect is not yet included. 

On the other hand Fig. shows our most interesting result, that is Higgs production 
cross section in the electron-positron collider including both QED ISR effect according to 
eq. (3.3) for ISR type (c) and for several examples of the machine energy spread parameter 
5 = 4.2,8,15,30,100MeV. 

In Fig. I^we plot again the same Higgs production cross section in the electron-positron 
collider including both QED ISR effect and machine energy spread in the narrower energy 
range and for the energy spread parameter 6 = 4.2MeV which is ambitiously aimed at the 
FCCee collider, and also for more realistic 6 = 8.0MeV. The ratio of cross section with 
respect to Born is also shown in the lower plot there. 

The same results as in Fig. are also shown in Tab. [^for three values of the energy, at 
the resonance peak E = Mh and near the peak E = Mh ± /%-, with the 4-digit precision. 
As we see there, the combined reduction of the Higgs production process at the resonance 
peak due to QED effect is 0.347 and goes down to 0.170 for the machine energy spread 
equal to Higgs width ~ 4.2MeV. 

The suppression factor of the peak cross section due to QED ISR can be also quite 
well reproduced, and in this way cross checked, with very simple approximate calculation 
presented in Appendix B. This approximation works even better in case of the weaker 
ISR effect for muon beams. 

The above results can be used as an input for further studies of the practical observ¬ 
ability of the Higgs resonant cross section in the future electron-positron colliders. 

In the above numerical exercises we have used discrete values of the machine energy 
spread S. In Fig. |^we also show the continuous dependence on 6 of the Higgs production 
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(7b 

(1) 

(2) 

(3) 

(11 

CTR 

O-R 

M 

O-R 

Eo 

67.8198 

37.14 

17.34 

11.06 

0.548 

0.256 

0.163 

Eo + Eb 

13.5636 

10.36 

13.62 

10.32 

0.764 

1.004 

0.761 

Eq — Eb 

13.5643 

7.40 

11.95 

9.62 

0.546 

0.881 

0.709 


Table 3: Cross section of the direct Higgs production process in the collider including 

effect of ISR type (c) and the machine energy spread. The same values of the energy spread 
6 an energy E are used as in Tabj^ The Monte Carlo integration error is below O.OSpb. 


process at the resonance peak, E = Mh, divided by the Born cross section, both for QED 
ISR switched on and off. More precisely, we are plotting there the following two ratios 


aBiMn) '''' asiMn) 


(4.1) 


for ISR off and on respectively. 

One could easily expand the above numerical exercises to more values in the 2- 
dimensional space of the energy E and machine energy spread 6, but we think that 
numerical results collected above are complete enough and can serve as a reliable starting 
point for all studies of the Higgs resonance observability in the future electron-positron 
colliders. 

It is worth to pointing out that because of the normalization of convolution as a 
function of 6 ai E = Mb to unity we can construct the following quantity that is close 
unity when 5/Tff is not too big. In Fig. [^the following quantity 


Rv{E,6) 


aniE) [a^^{E,0)) 


(4.2) 


is plotted. If Rv{E,S) ~ 1 then double convolution can be replaced by the following 
handy approximation 

_CO 

^conv / TT' r\ ^conv / rp n\ B 

^(c) [E, 0) ~ (F, 0) — 

As seen in Fig. [^the above approximation for E = Mb works reasonably, within 10% for 
6 < Eb, but its validity deteriorates signihcantly for S > Fb and in such a case numerical 
double convolution is unavoidable. 



5 Numerical results for muon colliders 

The calculations of the previous section can be easily extended to the case of muon collider, 
in the so-called Higgs factories (for overview see, e.g., [^), by means of replacing irie with 
rrifj, and B^e with in the equations of section]^ (see also or 0)' This kind of colliders 
will produce clean sample of Higgs boson without much background and therefore would 
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allow to measure Higgs mass and its properties very precisely. Our analysis is vital in view 
of MAP (Muon Accelerator Program) at Fermilab and the new ideas of producing 
muon monochromatic beams, that would allow to limit the major obstacle, the machine 
energy spread. 

Thanks to much higher branching ratio for the Higgs to muon pair, the production 
cross section for the H is much higher. Also, QED effects for muon collider are 

roughly factor 2 weaker, simply because for muon beams at the Higgs peek 7 = 0.0611, 
as compared to 7 = 0.1106 for electron beams. 

Numerical results for muon collider are presented in Fig. and Tab. and correspond 
to Higgs production cross section dependence on energy E and machine energy spread S 
shown in Fig. |^and in Tab. As seen in this table, reduction factor of the cross section 
at the resonance peak due to QED ISR is now only 0.548. For the machine energy spread 
6 = 8 MeV it deteriorates down to 0.162 (which is accidentally comparable to electron 
collider case with 5 = 4MeV) and for 6 = 4.2MeV it is equal 0.27^ 

6 Summary 

It was analysed numerically what is the influence of the QED ISR and machine energy 
spread on the the resonant Higgs boson production cross section, the so called Higgs 
line-shape, in the process e+e“ —)■ H and —)■ H. 

It was found that for electron collider the QED ISR reduces by itself peak cross section 
by factor 0.338. The QED higher order uncertainty of this results was estimated to be 
below 0.3%. The proper double convolution of the QED radiative spectrum with the 
machine energy spread was performed and doubly cross-checked. For instance machine 
energy spread of the same size as Higgs width (6 ~ 4MeV) reduces Higgs production cross 
section further down to 0.172 cts. These results are compatible (albeit slightly different) 
from these in the preliminary analysis shown at the recent FCCee workshop, see [^. It 
forms a solid basis for any analysis of the observability of the Higgs resonant cross section 
in the future e’''e“ colliders. 

The same analysis was also repeated for the muon colliders, where the reduction factor 
due to QED ISR at the peak position was found to be 0.548 and, for instance, it gets 
reduced further down to 0.162 for the energy beam spread being twice the Higgs width. 


®This agrees only rather crudely with 0.25 result in ref. for (5 
analytical formula. 


4MeV, obtained from approximate 
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A Non-uniqueness of the Breit-Wigner formula 

In this short Appendix we want to comment on the issue of the form of the non-uniqueness 
of the Breit-Wigner formula (2.1) and we shall show that for extremely narrow Higgs boson 
resonance this non-uniqueness is numerically completely irrelevant. 

The literature on this question is numerous, for instance in the context of the preci¬ 
sion measurements of the Z and W boson production and decay the reader may consult 
refs. 21 24 . In particular the use of the complex energy poles of the propagator of virtual 


particles as a natural solution was advocated, see ref. 24 


Following ref. (eqs. (1.6) through (1.9)), the Born term corresponding to (2.1) is 

. ^ s 


o-b(s) 


(s-M^)2 + G'2' 


(A.l) 


In case of constant width G = Mjj'jH = 527.94 where 'Jh = ^ = 3.34 ■ 10 while for 
the s dependent one we have G = G{s) = s^h- The maximum of the cross section is at 

(A.2) 


^ = Eo = Mh{1 + ^ Mh{1 + '^) = 125.700000035 GeV 


for the hxed width and 


2 

Eo = Mh{1 + ~ Mh( 1 - ^) = 125.69999996 GeV 


(A.3) 


for s-dependent width. The difference between Eq in the above two cases may be regarded 
as 7 X 10“® GeV. To our knowledge it is practically impossible to achieve energy resolution 
in accelerators of that order of magnitude. 

Finally, the formula with the so-called running width (for different variants of Breit- 
Wigner see, e.g., 25 and references therein) reads 

1 


o-B(sj 


(s-M^)2 + G2’ 
where G = s'Jh- The maximum of the cross section is at 

Mh 


■y/Tfi- — 








(A.4) 


(A.5) 


which, as above, is beyond measurability in any future accelerator experiment. 


B Approximate formulas for QED ISR 

Starting from eq. ( |2.2[ ) we are going to provide a very simple approximate formula that 
allows to estimate the suppression factor of the peak cross section at s = Mfj due to QED 
ISR corrections, thus providing quick/easy cross-check of more sophisticated numerical 
calculations. 














In the most simplified soft photon approximation ISR formula of eq. (2.2) aX E = Mjj 
reads 




H 


'yv'^ — v))dv. 


(B.l) 


It is well known fact that Breit-Wigner profile drop sharply around |s — M^| = R 


H 


and one may therefore approximate it by the following rectangular shape 

(Tb{.s) ~ cts(M^)0(|s - M|| < r|), 

which translates into v < Eh/Mh integration limit. As a result we obtain the following 
approximate suppression factor for ISR corrections 


risR = 


’’liMl) 1 f 


'yv'^ ^(TB{M‘fj)dv = 


r, 


H 


M, 


H 


(B.2) 


For e^e colliders with 7 = 0.1106 the ISR suppression factor is then estimated to be 


nsR = 0.3199. (B.3) 

This agrees reasonably well with the exact ISR factor of 0.347, seen in Tab. The 
relative error of the approximate formula is about 8% For collider with smaller 

7 = 0.0611 approximate ISR suppression factor 


rjsR = 0.5329, (B.4) 

is obtained, to be compared with the corresponding value of 0.548 in Tab. The relative 
error of the approximate formula is now merely 3%. 
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Figure 1; Study of the pure QED effect in the Higgs line-shape for electron collider. The 
plots show Born cross sections of eq. (2.1) ai|i^ cross section affected by QED ISR, following 
eq. (2.2) for three types I = a, b,c of the QED radiator functions defined in eq. (2.3). The 
ratios with respect to Born cross section are also shown. 

















Figure 2: Higgs production cross section in the electron-positron collider for several values of 
the machine energy spread ^ = 0,4.2, 8,15,30, lOOMeV. The QED ISR effect is not included. 
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Figure 3: Higgs production cross section in the electron-positron collider for several values 
of the machine energy spread S = 0,4.2, 8,15, 30, lOOMeV. The QED ISR effect is included 
according to eq. (3.3), for ISR type (c). 
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Figure 4: Line-shape of the direct Higgs production process in the electro-positron collider 
including effect of ISR type (c) and the machine energy spread. No energy spread is included 
(only ISR is on) in the line marked with (1), while lines marked with (2) and (3) include ISR 
and energy spread 6 = 4.2 MeV and 6 = 8.0 MeV according to eq. (3.3). Reference Born 
cross section is also shown and the ratios with respect to Born are also plotted. 
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Figure 5: The dependence of the Higgs production process at the resonance peak, E = Mh 
in the electron-positron collider on the energy spread 5, divided by the Born cross section, for 
QED ISR switched on and off. 
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Figure 6: Line-shape of the direct Higgs production process in the collider including 

effect of ISR type and the machine energy spread. The same input and notation as in Fig. 
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Figure 7: Numerical cross-check of the validity of the approximation shown in eq. (4.2) for 


E = Mh 
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